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Abstract. We present deep GMRT 21cm absorption spectra of 10 damped Lyman-a systems (DBAs), of which 8 are 
at redshifts z ^ 1.3. HI absorption was detected in only one DLA, the a = 0.5318 absorber toward PKS 1629+12. 
This absorber has been identihed with a luminous spiral galaxy; the spin temperature limit (Tg < 310 K) derived 
from our observations continues the trend of DBAs associated with bright spirals having low spin temperatures. 
In seven of the remaining 9 systems, the observations place strong lower limits on the spin temperature of the HI 
gas. 

We combine this sample with data taken from the literature to study the properties of all known DBAs with 
21cm absorption studies. The sample of DBAs which have been searched for 21cm absorption now consists of 31 
systems, with Ts estimates available in 24 cases; of these, 16 are at 2 < 2 and 8 at 2 > 2, with 11 (all at 2 < 1) 
having optical identifications. For the latter 11 DBAs, we find that all of the low Ts DBAs have been identified 
with large, luminous galaxies, while all the DBAs with high spin temperature {Ts ^ 1000 K) have been identified 
either with BSBs or dwarfs. Further, we find no correlation between impact parameter and spin temperature; it 
is thus unlikely that the high measured Ts values for DBAs arise from lines of sight passing through the outskirts 
of large disk galaxies. Instead, the spin temperature of DBAs appears to correlate with the host galaxy type. 
The trend (noted earlier by Ghengalur & Kanekar 2000) that low 2 DBAs exhibit both high and low Ts values 
while high redshift (2 ^ 3) DBAs only show high spin temperatures is present in this expanded data set. Based 
on this difference in spin temperatures, the Gehan test rules out the hypothesis that DBAs at 2 > 2 and DBAs 
at 2 < 2 are drawn from the same parent population at ~ 99% confidence level. 

Finally, we use the new GMRT spectra along with 2 spectra from the literature to estimate upper limits on the 
fraction of cold HI, /cnm, in DBAs at 2 ^ 3. For local spirals, /cnm ~ 0.5; in contrast, we find that /cnm < 0.3 
in all 7 high 2 absorbers, and /cnm < 0.1 in 5 of the 7 cases. 
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1. Introduction 

While molecular gas can be detected in emission at high 
redshifts (albeit only in a handful of somewhat unusual 
systems), it is as yet a severe challenge to detect neutral 
atomic gas in emission, even at redshifts as low as 2 ~ 0.2. 
As a result, absorption spectra against bright background 
sources provide the sole source of information on neutral 
atomic gas in the high redshift Universe. An important 
advantage of using absorption spectra as a probe of sys¬ 
tems lying along the lines of sight to distant QSOs is that 
such studies tend to be unbiased, at least to first order - 
sources which give rise to absorption are presumably typ- 
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ical systems at that redshift. In contrast, emission studies 
of a flux-limited, high redshift sample are usually domi¬ 
nated by the brightest (often atypical) sources. 

The large absorption cross-section of the Lyman- 
a transition and the dominance of hydrogen in baryonic 
material makes the Lyman-a line one of the easiest to de¬ 
tect in QSO spectra. In fact, lines arising from neutral gas 
with column densities as low as 10^^ atoms cm“^ can be 
routinely detected with present day telescopes. However, 
although these lines (the Lyman-a “forest”) are by far the 
most numerous, they contain only a small fraction of the 
total neutral gas at high redshifts; most of the neutral gas 
is, in fact, contained in relatively rare, high column density 
(.^Hi ^ 10^° atoms cm“^) systems, the so-called damped 
Lyman-a absorbers (DLAs). As the major repository of 
neutral gas at high redshift, these systems are natural can- 
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didates for the precursors of today’s galaxies. However, de¬ 
spite their importance in the context of galaxy evolution, 
the nature of high redshift DLAs is presently quite con¬ 
troversial, with models ranging from large rapidly rotat¬ 
ing disks ( Wolfe et al. 1986 ; Prochaska & Wolfe 1997) to 
merging sub-galactic blobs ( Haehnelt et al. 1998D to out- 
flows from dwarf galaxies (Schaye 2001) present in the re¬ 
cent literature. 

A successful model for DLAs needs to not only reproduce 
the observed Hg( 2 ) and column density distribution but 
also to describe the observed physical conditions in the ab¬ 
sorbing gas. The latter is important for our understanding 
of DLAs, quite independent of their precise size and struc¬ 
ture. Conversely, knowledge of these conditions is useful in 
modelling the further temporal evolution of damped sys¬ 
tems. HI 21cm observations (possible only for those DLAs 
for which the background quasar is also radio loud) are 
particularly useful because they allow one to estimate the 
average spin temperature Tg of the absorbing gas. DLA 
spin temperatures have been, in general, found to be sig¬ 
nificantly higher than those typically seen in lines of sight 
through the Galactic disk, or in those through the disks 
of nearby galaxies (see e.g . Chengalur fc Kanekar 200C , 
Kanekar & Chengalur 2001 and references therein). The 
only exceptions to the above are DLAs which have been 
identified as associated with low redshift (2 < 0.6) 

spiral galaxies - these systems all have low spin tem¬ 
peratures, comparable to those seen in the Milky Way 
and nearby spirals. Unfortunately, the number of sys¬ 
tems for which deep 21cm observations are available 
is relatively small; further, while the present data 
( [Chengalur fc Kanekar 200'o| ) suggest that high 2 DLAs 
are different from low 2 ones, the observations are cur¬ 
rently biased toward systems at 2 < 1, with not too many 
systems of the sample at high redshift. Are high redshift 
DLAs indeed systematically different from low redshift 
ones ? Also, how signihcant is the fact that only DLAs as¬ 
sociated with spiral galaxies have low spin temperatures ? 
To find answers to these and similar questions, it is crucial 
to increase the number of systems with 21cm observations, 
particularly at high redshifts. We present, in this paper, 
deep Giant Meterwave Radio Telescope (GMRT) observa¬ 
tions of 10 DLAs, 8 of which are at 2 > 1. Two of these 
systems, the 2 = 2.9084 DLA towards TXS 2342-1-342 and 
the 2 = 3.0619 absorber towards PKS 0336—014 have been 
earlier observed by Carilli et al. (1996); the present obser¬ 
vations are far more sensitive than the earlier ones. The 
observations themselves are summarised in Sect. while 
the results are presented in Sect. ^ and their implications 
for physical conditions in DLAs discussed in Sect. ||. 


2. Observations and Data Analysis 

2.1. The GMRT observations 

The GMRT has hve currently operational frequency 
bands, at 150 MHz, 233 MHz, 327 MHz, 610 MHz and 
1420 MHz. The present observations utilised the 327, 


610 and 1420 MHz bands, which cover redshift ranges 
2 ~ 2.9 — 3.6, 2 1.15 — 1.5 and 2 ~ 0 — 0.6 respec¬ 

tively, with reasonable sensitivity (i.e. better than half 
the sensitivity at the band centre). The observations were 
carried out between January and November 2001, using 
the 30-station FX correlator as the backend; this gives a 
fixed number of 128 channels over a bandwidth which can 
be varied between 64 kHz and 16 MHz. All our observa¬ 
tions (except for those of PKS 2128—123 (bandwidth = 
2 MHz) and PKS 1629-1-12, which are described below) 
used a bandwidth of 1 MHz, sub-divided into 128 chan¬ 
nels; this gave a channel resolution of 7.8 kHz. The num¬ 
ber of available antennas varied between 18 and 26, due 
to various debugging and maintenance activities. One (or 
more) of the standard calibrators 3G48, 3C147 or 3C286 
was used in all cases to calibrate the absolute flux scale 
and the system bandpass; these sources were observed at 
least every two hours in all cases, to ensure a good band¬ 
pass calibration. A nearby standard phase calibrator (see 
Gols. 4 and 5 of Table 1) was observed every forty minutes 
in all cases to determine the instrumental phase, except in 
the case of PKS 2128—123 which is itself a phase calibra¬ 
tor for the GMRT. Our experience with the GMRT indi¬ 
cates that the flux calibration is reliable to ~ 15%, in this 
observing mode. Observational details are summarised in 
Table 1, where the sources are arranged in order of increas¬ 
ing redshift; note that the values listed here for spectral 
resolution and RMS noise are before any smoothing. 

The only source towards which we detected 21cm ab¬ 
sorption, viz. PKS 1629-1-12, was initially observed with 
a 1 MHz bandwidth on the 12th and 22nd of October 
2001, with 20 antennas and a total on-source time of 
~ 6 hours. These observations resulted in the detection of 
a narrow 21cm absorption feature; the source was hence 
re-observed on the 17th and 18th of November, 2001 with 
a 0.25 MHz bandwidth (yielding a resolution of 1.95 kHz, 
i.e. 0.6 km s“^), to try to resolve out the absorption line. 
The total on-source time of these latter observations was 
13 hours, with 25 antennas. In all cases, 3G286 was used 
for absolute flux and system bandpass calibration while 
the compact source 1640-1-123 was used for phase calibra¬ 
tion. 


2.2. Data Analysis 

The data were analysed in classic AIPS using standard 
procedures. After the initial flagging and gain and band¬ 
pass calibration, continuum images were made for all 
sources and then used to self-calibrate the U-V data. This 
was carried out in an iterative manner until the qual¬ 
ity of the image was found to not improve on further 
self-calibration. In general, a few rounds of phase self¬ 
calibration were carried out on all fields, followed by one or 
two rounds of amplitude self-calibration. The continuum 
emission was then subtracted from the multi-channel U- 
V data set, using the AIPS task UVSUB; any residual 
continuum was subtracted by fitting a linear baseline to 




















Kanekar & Chengalur: 21cm absorption in high z DLAs 


3 


Table 1. Observing details 


QSO 

'^•abs 

i^obs 

MHz 

Phase 

Cal. 

Phase Cal. 
Flux 

Jy 

No. of 

antennas 

Time 

Hrs. 

Source 

Flux 

Jy 

Channel 
resolution 
km s“^ 

RMS 

mJy 

PKS 2128-123 

0.4298 

993.43 

a 

a 

20 

8 

1.90 

4.7 

1.4 

PKS 1629+12 

0.5318 

927.28 

1640+123 

2.5 

25 

13 

2.35 

0.6 

2.8 

PKS 0215+015 

1.3439 

605.72 

0204+152 

5.0 

26 

5 

0.92 

3.9 

2.2 

QSO 0957+561A 

1.3911 

594.04 

0834+555 

7.4 

26 

6.5 

0.59 

3.9 

2.4 

PKS 1354+258 

1.4205 

586.82 

3C286 

b 

22 

5.5 

O.SO'^ 

4.0 

2.8 

TXS 2342+342 

2.9084 

363.42 

3C19 

8.3 

23 

8 

0.31 

6.4 

2.8 

PKS 1354-107 

2.966 

358.15 

3C283 

21 

22 

9.5 

0.12 

6.5 

2.0 

PKS 0537-286 

2.974 

357.43 

0521-207 

7.2 

25 

5 

1.05 

6.5 

2.4 

PKS 0336-014 

3.0619 

349.69 

0323+055 

7.4 

26 

5 

0.94 

6.7 

2.1 

PKS 0335-122 

3.178 

339.97 

0409-179 

6.1 

18 

5.5 

0.68 

6.9 

1.9 


“ No phase calibrator was used, as PKS 2128—123 is itself a phase calibrator for the GMRT. 

3C286 was also used as the flux calibrator; we hence do not have an independent estimate of its flux. 
PKS 1354+258 is resolved by the GMRT synthesised beam; this is the peak flux. 


the U-V visibilities, using the task UVLIN. The contin¬ 
uum subtracted data were then mapped in all channels 
and a spectrum extracted at the quasar location from 
the resulting three dimensional data cube. In the case 
of PKS 1629+12, which showed 21cm absorption, spectra 
were also extracted from other locations in the cube to en¬ 
sure that the data were not corrupted by RFI. In general, 
all the observing frequencies were found to be reasonably 
free of RFI, except for sporadic interference on individual 
antennas which was edited out. In the case of multi-epoch 
observations of a single source, the spectra were shifted 
to the heliocentric frame outside AIPS and then averaged 
together. 

3. Results 

The final absorption spectra towards the ten sources of our 
sample are shown in Fig. |^, arranged in order of increasing 
redshift; the measured source fluxes, channel resolutions 
(in km s“^) and RMS noise values are listed in Cols. 8, 9 
and 10 of Table 1. 

Our only detection of HI absorption is in the z = 
0.5318 DLA towards PKS 1629+12. Fig. ^ shows a more 
detailed version of the GMRT 0.25 MHz HI spectrum 
towards this source; the spectrum has a resolution of 


^ 0.6 km s“^ and an RMS noise of 2.8 mJy. Highly 
asymmetric HI absorption can be clearly seen, extending 
over a velocity width of ^ 40 km s“^. The peak absorp¬ 
tion occurs at a heliocentric frequency of 927.287 MHz, 
i.e. at a heliocentric redshift of 0.531787 ± 0.000002. 
This lies in between (and in reasonable agreement with) 
the r edshifts of the da mped Lyman-a line (zdla = 
0.532; Nestor et al. 2001) and the Mgl absorption (zMgi = 
0.5316; ( Aldcroft et al. 1994 ; Barthel et al. 1990|) . (Note 
that Aldcroft et al. (1994) quote z = 0.5313 for the strong 
Mgll and Fell absorption lines.) The flux of PKS 1629+12 
was measured to be 2.35 Jy; the peak optical depth is 
thus Tmax = 0.039. The equivalent width of the profile is 
f rdV = 0.494 ± 0.001 km s“^, obtained by integrating 
the observed absorption profile. 


4. The spin temperature 


In the case of an optically thin, homogenous cloud in ther¬ 
mal equilibrium, its HI column density, IVhi, 21cm optical 
depth, T, and spin temper ature, Tg, are related by the ex¬ 
pression (e.g. Rohlfs lO"^) 


A'hi = 


1.823 X IQi^Tg 

7 


rdP , 


( 1 ) 
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Fig. 1. GMRT HI absorption spectra towards the ten sources of our sample, arranged in order of increasing redshift. 
[A] 2 = 0.4298 DLA towards PKS 2128-123. [B] 2 = 0.5318 DLA towards PKS 1629+12 [C] 2 = 1.3449 DLA 
towards PKS 0215+015. [D] 2 = 1.3911 DLA towards QSO 0957+561A. [E] 2 = 1.4205 DLA towards PKS 1354+258. 
[F] 2 = 2.9084 DLA towards TXS 2342+342. [G] 2 = 2.966 DLA towards PKS 1354-107. [H] 2 = 2.974 DLA towards 
PKS 0537-286 [I] 2 = 3.0619 DLA towards PKS 0336-014. [J] 2 = 3.178 DLA towards PKS 0335-122. 


where Ahi is in cm“^, Tg in K and dV in km s“^. / is the 
covering factor of the absorber. Of course, a given line of 
sight through an absorber is quite likely to pass through 
a number of clouds, with different spin temperatures and 
column densities. In such a situation, the spin temperature 
obtained using the above expression is the column density 
weighted harmonic mean of the spin temperatures of the 


individual clouds (provided, of course, that all the clouds 
are optically thin). 

It should be noted that Eq. (|^ relates the HI col¬ 
umn density to the integrated optical depth J rdP; this 
means that, in the case of a non-detection, one needs to 
know the shape of the absorption profile to compute a 
lower limit on the spin temperature. This shape is, of 
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Table 2. The spin temperature for the 10 DLAs observed here 


QSO 

^ahs 

TVhi X 10^° 
cm“^ 

AF “ 
km 

RMS= 

mjy 

d 

"^max 

(l//)/rdF^ 
km s ^ 

Tg 

K 

Refs.-^ 

PKS 2128-123 

0.4298 

0.25 + 0.06 

7 

1.2 

< 0.0019 

< 0.014 

> 980 

1 

PKS 1629+12 

0.5318 

2.8 

0.6 

2.8 

0.039 

0.494 + 0.001 

310 

2 

PKS 0215+015 

1.3439 

0.8 

7.8 

1.6 

< 0.0052 

< 0.043 

> 1020 

3 

QSO 0957+561A*’ 

1.3911 

2.1 ±0.5 

- 

- 

- 

- 

- 

4 

PKS 1354+258*’ 

1.4205 

32 + 2 

- 

- 

- 

- 

- 

4 

TXS 2342+342 

2.9084 

20 + 0.5 

15 

2.0 

< 0.0192 

< 0.306 

> 3585 

5 

PKS 1354-107 

2.966 

6.0 

6.5 

2.0 

< 0.05 

< 0.345 

> 955 

6 

PKS 0537-286 

2.974 

2.0 

10 

1.9 

< 0.0054 

< 0.058 

> 1890 

6 

PKS 0336-014 

3.0619 

16 + 1.3 

20 

1.4 

< 0.0045 

< 0.095 

> 9240 

7 

PKS 0335-122 

3.178 

6 + 1 

15 

1.2 

< 0.0053 

< 0.084 

> 3920 

6 


“AL is the FWHM used to compute the spin temperature in the case of non-detections (see the text for more details). For 
1629+12, it is the velocity resolution of the observations. 

*'The large uncertainty in the covering factor makes it difficult to estimate the spin temperature (see the text for more details). 
■^This is the RMS acually measured from the smoothed spectra. The spectra were smoothed to a velocity resolution given in 
the preceeding column, except for the following cases for which the RMS was measured at a slightly better spectral resolution 
than AF, viz 12.8 km s“^ for TXS 2342+342, 13.4 km s"^ for PKS 0336-014 and 13.8 km s"^ for PKS 0335-122. 

‘^All limits (in this and other columns) correspond to 3 ct. 

®The covering factor / has been set to unity in all cases (see the text for d etails). _ 

■^ References for the H f column densitie s : 1. [Ledoux et al. 20 03; 2. iNestor et al. 200l|; 3. 

White et al. 1993; 6. 


4. Rao & Turnshek 200C; 5. 


Ellison et al. 2001a; 7. Prochaska et al. 2001 


Lanzetta et al. 1995 


course, not a priori known. However, since our limits on 
the optical depth imply that the HI line is optically thin, 
the line can be assumed to have a Gaussian shape; here, 
f rdV = l.ObTmaxAF, where r^ax is the peak optical 
depth (or the limit on the peak) and AV is the FWHM of 
the Gaussian. Further, for a gas cloud in thermal equilib¬ 
rium, the kinetic temperature Tk (with Tg « Tk, for a sin¬ 
gle cloud) and the FWHM of the observed line are related 
(Tk ~ 21.855AF^, with Tk in K and AF in km s“^); care 
must hence be taken to derive the limits on Tg in a self- 
consistent manner, so that the velocity width over which 
the integral is carried out is consistent with (at least) ther¬ 
mal broadening of the gas at the derived temperature (any 
bulk kinematic motion will further broaden the absorption 
features; we will assume that no such motions are present). 
For example, if the Scr lower limit on the spin temperature 
is 10000 K, one should use AF ^ 20 km s“^, while, if the 
limit is 1000 K, it is appropriate to use AF ^ 7 km s“^. 
We have, in the following analysis, used AF « 7, 15 and 
20 km s“^ for temperatures Tg ~ 1000, 5000 and 10000 K, 


respectively, and have smoothed the spectra to the above 
resolutions while computing the spin temperature. 

The HI column density of a DLA can be directly es¬ 
timated from the equivalent width of the Lyman-a pro¬ 
file. The original criterion for a system to be classified as 
a damped absorber (iVni > 2 x 10^° cm“^) was an ob¬ 
servational one, linked to the ability to identify damped 
profiles in medium resolution spectra from the relatively 
small telescopes used to carry out the early DLA surveys. 
Throughout this paper, however, we will consider all sys¬ 
tems for which the Lyman-a profile shows damping wings 
as DLAs, regardless of whether or not their HI column 
density is greater than the traditionally used threshold. 
We hence retain the absorbers towards PKS 0215+015 
and PKS 2128—123 in our DLA sample, despite their HI 
column densities being less than 10^° cm“^. 

Given the HI column density, a search for 21cm ab¬ 
sorption directly yields the spin temperature of the ab¬ 
sorbing gas (or limits on the temperature, in the case of 
non-detections of absorption), provided the covering fac- 
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Fig. 2. High resolution (0.6 km s“^) GMRT HI absorp¬ 
tion spectrum towards PKS 1629-1-12. Asymmetric HI ab¬ 
sorption is clearly visible, with the peak optical depth at 
2 = 0.531787. 

tor / is known. Unfortunately, the radio continuum from 
quasars tends to arise from a far more extended region 
than the optical or UV emission and it may thus transpire 
that some fraction of the radio continuum may “leak out” 
from around the cloud and result in an incorrectly esti¬ 
mated spin temperature. Similarly, the line of sight along 
which the 21cm optical depth has been measured need 
not be the same as that for which one has an estimate for 
the HI column density from the Lyman-a line. The frac¬ 
tion of radio emission originating in compact components 
spatially coincident with the UV point source can usually, 
however, be measured from VLBI observations (when such 
observations are available at sufficiently low frequencies). 
One can then estimate / and hence, the spin temperature; 
this is done next, for the ten sources of this paper, in order 
of increasing right ascension. The results are summarised 
in Table 2 in which Cols. 4 to 8 contain, respectively, the 
velocity width AV used to compute Tg, the spectral RMS 
over this (or a slightly smaller) velocity width, the peak 
optical depth Tmax, (1//) / ^dU, and the spin temperature 
Ts (or limits on the last three quantities). 

1. VLBA maps of PKS 0215-b015 at 2.3 and 8.4 
GHz (from the Radio Reference Frame Image Database 
(RRFID)) show that the source is exceedingly compact, 
with most of the flux contained within the central 30 milli- 
arcseconds. The quasar also has a flat spectrum between 
8.4 GHz and 610 MHz, again implying that it is quite 
likely to be compact even at the lower frequencies. It 
is hence likely that the covering factor / is close to 
unity. The HI colum n density of the z = 1.345 DLA is 
.^Hi = 8 X 10^® cm“^ ( Lanzetta et al. 1995 ); the spin tem¬ 


perature of the absorber is then Ts{3a) > 1020 K (using 
AU = 7.8 km s'^). 

2. The 2 = 3.178 DLA towards PKS 0335-122 is 
the highest redshift system of our sample, with Vhi = 
6 ± 1 X 10^° cm“^ (Ellison et al. 2001a). Unfortunately, 
no VLBI information exists in the literature on the back¬ 
ground quasar; however, this source has a fairly flat spec¬ 
trum between 5 GHz and 1.4 GHz, with fluxes of 380 mJy 
and 420 mJy at these two frequencies. This spectral flat¬ 
ness implies that the flux is likely to arise from a fairly 
compact source component, which is hence likely to be 
covered by the DLA. This gives a lower limit to the cov¬ 
ering factor, / > 0.7, which yields Ts(3ct) > 2740 K; if 
/ = 1, one has Ts(3ct) > 3920 K (AU = 15 km s“^). 

3. In the case of PKS 0336—014, 5 GHz VLBI observa¬ 
tions ( iGurvits et al. 19^ ) have shown that the source is 
highly compact at this frequency, with the entire flux con¬ 
tained within the central ~ 7 milli-arcseconds. However, 
the source also has a negative spectral index, indicat¬ 
ing the presence of some extended emission. In the ab¬ 
sence of VLBI information at lower frequencies, we will 
assume a covering factor of unity. The HI column den¬ 
sity of the 2 = 3.0619 DLA is IVhi = 1.6 ± 0.13 x 
10^^ cm“^ (Prochaska et al. 2001) giving a spin temper¬ 


ature Ts(3ct) > 9240 K (AU = 20 km s“^). Carilli et al. 
(1996) estimated Tg > 2500 K for this DLA. 

4. PKS 0537—286 has an exceedingly flat spectrum 
between 8.4 GHz and 327 MHz; for example, the QSO 
flux density was measured to be 0.9 Jy at 8 GHz (Parkes 
survey, [Wright & Otrupcek 1990[), 1.2 Jy at 4.85 GHz 
( priffith et al. 1994D, 838 mJy at 1.4 GHz (NRAO-VLA 
Sky Survey, Condon et al. 1998) and 1.05 Jy in our GMRT 
327 MHz image. This implies that the source is likely 
t o be highly compac t. VLBA observations at 5 GHz 
( Fomalont et al. 200C ) have shown that the flux within 
the central 10 milli-arcseconds is ~ 500 mJy, imply¬ 
ing a covering factor greater than 0.5. In fact, these 
observations recovered the entire 5 GHz flux of 1.2 Jy 
on the shortest VLBA baseline (2 x 10® wavelengths, 
Hirabayashi et al. 200(]| ), implying that the entire source 
is contained within ~ 0.1 arcseconds. It is thus quite 
likely that the covering factor is close to unity; we will 
use / = 1 in this analysis. The HI column density 
of the 2 = 2.9742 DLA has been found to be 2.0 x 
10^® cm“^ (Ellison et al. 2001a); we then obtain a 3a 


lower limit of 1890 K on the spin temperature of the ab¬ 
sorbing gas (AU = 10 km s“^). 

5. VLBI mapping of QSO 0957-1-561 has been carried 
out on a number of occasions, at various frequencies at 


and above 1.6 GHz (see, for example, Garrett et al. 1994, 
pampbell et al. 1995 ). The flux of component A has 
been measured to be only about 32 mJy at 1.6 GHz 
(Garrett et al. 1994); clearly, most of the flux seen in 
the 610 MHz GMRT map comes from extended emission 
which is unlikely to be covered by the DLA. The cover¬ 
ing factor is thus quite uncertain for this absorber, al¬ 
though it is likely to be low. Given this, limits on the spin 
temperature derived from the present data are unlikely to 
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be reliable. In fact, it is known that the DLA does not 
cover the second image (B) of this gravitational lens sys¬ 
tem ( Rao fc Turnshek 2000D . We will hence not include 
this system in the later analysis and also do not list limits 
on the spin temperature and the optical depth in Table 2. 

6. PKS 1354-1-258 is also an extended source at 

610 MHz, with no VLBI information available in the lit¬ 
erature; the DLA covering factor is hence again quite 
uncertain. However, the extremely high column den¬ 
sity of the absorber (A^hi = 3.2 ± 0.2 x 10^^ cm“^; 

Rao & Turnshek 2000) results in high lower limits on the 
spin temperature even in the case of low covering factors. 
For example, / = 0.1 gives Ts{3a) > 1100 K, while / = 1 
would imply Ts{3a) > 7000 K (assuming velocity widths 
of 8 km s“^ and 20 km s“^ respectively. We will however 
again not quote a spin temperature in Table 2 due to the 
uncertainties in the covering factor. 

7. No information exists in the literature regarding the 
VLBI structure of PKS 1354—107. However, the source 
has a strongly inverted spectrum and is thus likely to be 
very compact. We will hence assume a covering factor of 
unity. Combining this with the HI column density of 6 x 
10^° cm“^ (Ellison et al. 2001a) gives Ts{3a) > 955 K 
(AV = 6.5 km s-^). 

8. The HI column density of the z = 0.5318 
DLA towards PKS 1629-1-12 is Vpi = 2.8 x 
10^° cm“^ (Nestor et al. 2001). PKS 1629-1-12 is a well- 


9. PKS 2128—123 has a very flat spectrum between 
8.4 GHz and 750 MHz, with most measurements between 
these frequencies yielding a flux density of around 2 Jy 
(see the NEDQ for details). VLBA maps at 2.3 and 8.4 
GHz show that most of this flux is contained within 
~ 20 milli-arcseconds, implying that the covering fac¬ 
tor is likely to be close to unity. The HI column den¬ 
sity of the z ~ 0.4298 absorber is Vhi = 2.5 ± 0.06 x 
10^® cm“^ ( Ledoux et al. 20^ ); the 3 (t lower limit on 
the spin temperature is then Ts{3(t) > 980 K (AV = 
7 km s-i). 

10. Finally, TXS 2342-1-342 has a fairly flat spec¬ 


trum between 1.4 GHz and 8 GHz ( Wilkinson et al. 1998 
[White fc Becker 1992), indicating that the quasar is very 


compact at these frequencies. However, as Carilli et al. 
(1996) point out, the spectrum rises slowly down to the 
redshifted 21cm frequency of 356 MHz, indicating the 
presence of a more extended source component. Carilli 
et al. (1996) argue that the DLA is likely to at least 
cover the compact component and the covering factor 
is hence larger than 0.5; it is, however, unclear whether 
the extended flux is also covered. We will use / = 1 
in the analysis (following parilli et al. 1996 ); the HI col¬ 
umn density of Nm = 2 ± 0.5 x 10^^ cm“^ then yields 
T,(3cr) > 3585 K (AV = 15 km s"!). Carilli et al. (1996) 
estimated Tg > 1800 K for this DLA. 


known compact steep spectrum source with an extremely 
one-sided structure ( ^aikia et al. 199Ci| ). The 2.3 GHz 
VLBI map of Dallacasa et al. (1998) shows two com¬ 
ponents with fluxes of 110 mJy and 150 mJy. Only the 
110 mJy western component is seen in their 8.4 GHz 
map, with almost the same flux; they hence identify 
this as the flat spectrum core of the source. The sep¬ 
aration between the two components is 1.14 arcsec (i.e. 
~ 6.4 kpc at z = 0.5318, using Hq = 65 km s“^ Mpc“^ 
and an n = 1, A = 0 cosmology). On the other hand, 
the core flux density was measured to be 179 mJy at 
408 MHz from a ^ 0.6 x 0.5 arcsec resolution Merlin im¬ 
age ( ^aikia et al. 1990 ). Using the 2.3 GHz and 408 MHz 
flux densities to estimate the core spectral index, we ob¬ 
tain a fliDc density of ~ 142 mJy at the redshifted 21cm 
line frequency of ~ 927 MHz. If the DLA covers only the 
core, the spin temperature of the absorber is Tg ~ 20 K, 
a remarkably small value. However, the 21cm absorption 
is spread over ^ 40 km s“^ and is hence unlikely to arise 
from a single line of sight. Further, the 408 MHz Merlin 
map ( ^aikia et al. 199C ) shows that the entire flux at this 
frequency is contained within ~ 3 arcsecs (~ 16.8 kpc 
at z = 0.5318). It is thus possible that the DLA covers 
some of the extended flirx of the background source, be¬ 
sides the core. VLBI observations in the redshifted 21cm 
line would be necessary to detect the source components 
towards which the absorption occurs, and hence, the cov¬ 
ering factor. For now, we note that a covering factor of 
unity results in Tg = 310 AT; this is a definite upper limit 
on the spin temperature, which must then lie between 20 K 
and 310 K. 


5. Discussion 

Neutral gas in the galaxy is believed to exist in two sta¬ 
ble phases, a cold dense phase (the cold neutral medium, 
GNM), with typical temperatures of ~ 100 K, and a warm 
rarefied phase (the warm neutral medium, WNM), with 
typical temperatures of ^ 8000 K. It should be empha¬ 
sised that such a multi-phase medium exists only over a 
finite range of interstellar pressures (Wolfire et al. 1995); 
at higher pressures only the CNM phase exists, while at 
lower pressures only the WNM phase is found. Further, 
although the multi-phase structure results from a bal¬ 
ance between the heating and cooling rates, which them¬ 
selves depend on various quantities such as the radia¬ 
tion field, the metallicity, the dust content, etc., it has 
been shown that a two-phase pressure equilibrium is pos¬ 
sible over fairly wide ranges of the above quantities (and, 
of course, the interstellar pressure; Wolfire et al. 1995). 


On the observational front, the existence of neutral gas 
in a multi-phase medium has been established not just 


locity clouds (e.g. 

Braun & Burton 2000t), dwarf galax- 

ies (e.K- 

Young & Lo 1997 

) and DLAs (Lane et al. 2000 

Kanekar et al. 2001b 

). In the latter case of the two DLAs 


at z = 0.0912 and z = 0.2212 toward B0738-I-313, it has 
been explicitly shown that the observed high spin temper¬ 
atures are a consequence of 70 — 80 % of the gas being in 
the warm phase. In this regard, the multi-phase structure 
of at least these two systems is more like that of dwarf 
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Table 3. Damped Ly-a Systems with 21cm Observations 


Name 

^ahs 

Whi 

10^° cm"^ 

/ 

(l//)JrdU 
km s ^ 

AU2i“ 
km s~^ 

S 

K 

Absorber 

ID 

5“' 

kpc 

Refs. 

0738+313 (A) 

0.09123 

15 + 2 

1 

1.0 + 0.01 

40 

825 + 120 

Dwarf 

< 4 

1,2,3 

PKS 0439-433“ 

0.10097 

7 

1 

0.07 + 0.02 

10 

7 

Spiral ? 

8.0 

4,5 

0738+313 (B) 

0.22125 

7.9 + 1.4 

1 

0.488 + 0.004 

30 

890 + 165 

Dwarf 

20 

1,3,6 

PKS 0952+179 

0.23779 

21 + 2.5 

0.25 

0.56 + 0.02 

30 

2055 + 330 

LSB 

< 7 

7,8,9 

PKS 1413+135“* 

0.24671 

200* 

0.12 

59.43 

39 

185 

Spiral 

7 

10,11,12 

PKS 1127-145 

0.3127 

51 + 9 

1 

3.074 + 0.002 

120 

910 + 160 

LSB 

< 10 

8,9,13 

PKS 1229-021 

0.39498 

5.6+ 0.5 

0.5 

1.81+0.01 

110 

170 + 15 

Spiral 

11 

14,15,16 

QSO 0248+430 

0.3941 

7 

7 

2.99 

40 

7 

- 

- 

17 

PKS 2128-123 

0.4298 

0.25 + 0.06 

1 

< 0.014 

- 

> 980 

- 

- 

18,19 

3C196 

0.43670 

6.3 + 1.6 

7 

7 

250 

- 

Spiral 

9.2 

14,20,21 

PKS 1243-072 

0.4367 

7 

1 

0.38 

23 

7 

Spiral 

11.3 

17,22 

3C446“ 

0.4842 

6.3 

0.1 

< 0.22 

- 

> 1570 

- 

- 

13,23 

AO 0235+164 

0.524 

50 

1 

13 + 0.6 

125 

210 

Spiral 

< 11.2 

24,25,26 

B2 0827+243 

0.5247 

2.0+ 0.2 

0.67 

0.33 + 0.03 

50 

330 + 70 

Spiral 

34 

7,8,9 

PKS 1629+12 

0.5318 

2.8 

1 

0.494 + 0.001 

40 

310 

Spiral 

17 

9,19 

PKS 0118-272“ 

0.5579 

7 

0.5 

< 0.12 

- 

7 

- 

- 

7,27 

3C286 

0.69215 

20 + 1.6 

1 ? 

0.91+0.09 

18 

1205 + 240 

LSB 

19 

15,28,29 

PKS 0454+039 

0.8596 

5 + 0.16 

1 

< 0.32 

- 

>855 

Dwarf 

6.4 

15,30,31 

PKS 0215+015“ 

1.3449 

0.8 

1 

< 0.043 

- 

> 1020 

- 

- 

19,32 

QSO 0957+561A 

1.3911 

2.1 + 0.5 

7 

- 

- 

- 

- 

- 

8,19 

PKS 1354+258 

1.4205 

32 + 2 

- 

- 

- 

- 

- 

- 

8,19 

MC3 1331+170 

1.77636 

15 + 1.5 

0.8 

0.5+ 0.1 

25 

1645 + 615 

- 

- 

31,33,34 

PKS 1157+014 

1.944 

63 + 10 

0.55 

4 + 0.2 

60 

865 + 190 

- 

- 

31,35,36 

PKS 0458-02 

2.03945 

45 + 8 

1 

6.4+ 0.4 

30 

385 + 100 

- 

- 

34,37,47 

PKS 0528-2505 

2.8110 

22 + 1 

1 

< 1.70 

- 

> 710 

- 

- 

38,39,40 

TXS 2342+342 

2.9084 

20 + 0.5 

1 

< 0.306 

- 

> 3585 

- 

- 

19,41 

PKS 1354-107 

2.966 

6 

1 

< 0.345 

- 

> 955 

- 

- 

19,42 

PKS 0537-286 

2.9742 

2 

1 

< 0.058 

- 

> 1890 

- 

- 

19,42 

PKS 0336-014 

3.0619 

16 + 1.3 

1 

< 0.095 

- 

> 9240 

- 

- 

19,41,46 

PKS 0335-122 

3.178 

6 + 1 

1 

< 0.084 

- 

> 3920 

- 

- 

19,42 

PKS 0201+113 

3.3875 

18 + 3 

1 

< 0.30 

- 

> 3290 

- 

- 

41,44,45 


* T he 2 = 0.24671 abso rber towards PKS 1413+135 may be an “associated” absorption syste m, rather than an intervening 
one ( Perlman et al. 2002| ). Further, its HI column density is estimated from an X-ray spectrum (Stocke et al. 1992). 

“ The velocity width AT quoted here is the entire velocity range over which absorption is seen. 

** The spin temperature values (or 3cr limits) quoted here have been consistently (re-)computed from the published data using 
E<1- ( 0 )> and the latest available values for all parameters in the equation. The values quoted here and the errors have been 
rounded off to the nearest 5 K. The errors on Nm and hence, on Ts, have been quoted when available in the original references; 
in the case of asymmetric errors, the larger value has been quoted. In certain instances, the values may differ slightly from that 
quoted by the authors in the original reference. The only exception is PKS 0201+113, for which multiple and quite different 
values are quoted in the literature. We use the value from Kanekar & Chengalur (1997), modihed by the new HI column density 
of 1.8 ± 0.3 X 10^*^ cm“^, from Ellison et al. (2001b). 

“ Candidate damped Lyman-a systems, based on lUE spectra, metal lines or X-ray spectra. 

The impact parameter b has been calculated using an D = 1, A = 0 Universe, with Hq = 65 km s“^ Mpc“^. 
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galaxies (see, for example, Young fc Lo 1997 ) than large 
spiral disks. 

Table 3 contains a list of the 31 known candidate and 
confirmed DLAs for which 21cm observations are avail¬ 
able in the published literature (this is an update of the 
list in [Chengalur fc Kanekar 20001 ). Note that the Lyman- 
a line has so far not been observed for the DLAs at 
z = 0.24671 towards PKS 1413-1-135, z = 0.3941 towards 
QSO 0248-f430 and z = 0.43669 towards PKS 1243-072; 
however, all three systems show strong 21cm absorp¬ 
tion, implying (for reasonable values of Ts) HI column 
densities higher than 10^° cm~^ (ICarilli et al. 1992 ; 
Kanekar et al. 20'0^ ; Lane fc Briggs 2002| ) and have hence 
been included in the present sample. Further, the z = 
0.24671 absorber towards PKS 1413-1-135 may be an 
associated system ( Perlman et al. 2002 ); the HI column 
density listed i n Table 3 is estim ated from the deficit 
of soft X-rays ( ^tocke et al. 1992 ). No similar direct es¬ 
timate exists for the HI column densities of the sys¬ 
tems towards PKS 1243-072 and QSO 0248-f430. We 
also note that two candidate DLAs (on the basis of 
metal lines), at z = 0.55 79 towards PKS 0118—272 
(Kanekar & Chengalur 2001) and z = 0.101 towards 
PKS 0439—433 (Kanekar et al. 2001a) have, for complete¬ 
ness, been included in this table; however, they will not 
be used in the later analysis, as their HI column densi¬ 
ties (and hence, spin temperatures) are unknown. Cols. 3, 
4, 5, 6 and 7 of the table contain the HI column density 
(in units of 10^° cm“^), the adopted covering factor /, 
the “true” equivalent width ((1//) / rdP), the velocity 
spread of 21cm absorption (between nulls, in km s“^) and 
the spin temperature. Besides these. Cols. 8 and 9 list the 
nature of the objects identified as the DLA hosts and their 
impact parameters (5, in kpc) to the QSO line of sight. Of 
course, it is possible that the latter identifications are in 
error in some cases if, for example, the absorption arises 
in an object closer to the QSO but too faint to be iden¬ 
tified in optical/IR images. With this caveat in mind, we 
note that seven of the DLA hosts have been identified with 
large, luminous galaxies, with luminosities close to L*. 

Fig. ^ shows a plot of the measured spin temperatures 
as a function of redshift for the 24 DLAs of the complete 
sample which have reliable estimates of the spin temper¬ 
ature (note that this plot includes the z = 0.24671 ab¬ 
sorber towards PKS 1413-1-135 which may be an associ¬ 
ated system). The seven systems in Table 3 which are not 
included in Fig. ^ are the two 21cm absorbers towards 
PKS 1243-072 and QSO 0248-b430 for which the HI col¬ 
umn density is unknown (see above), the DLAs towards 
PKS 1354-1-258 and QSO 0957-1-561 for which the covering 
factor is highly uncertain (see Sect. the two candidate 
DLAs towards PKS 0118-272 and PKS 0439-433 and, 
finally, the z = 0.437 absorber towards 3C196, where the 
lines of sight to the optical QSO and the radio continuum 
trace very different paths through the absorbing galaxy. 
The dotted line in Fig. ^ is at 350 K - more than 80 % 
of spin temperature measurements in M 31 and the Milky 
Way have values < 350 K ( [Braun fc Walterbos 1992 ). It 



0 12 3 


Redshift, z 

Fig. 3. The spin temperature (LogjTg]) of the 24 DLAs 
of the complete sample with Tg estimates as a function 
of redshift z. Note that seven DLAs listed in Table 3 
have not been included here for reasons listed in the text. 
Objects identified with spiral galaxies are plotted using 
squares; the filled square shows the z = 0.24671 spiral 
towards PKS 1413-1-135, which may be an associated sys¬ 
tem. Other detections are shown as filled triangles, with 
non-detections shown as open triangles, with arrows. 


can be seen that the spin temperatures of the majority of 
DLAs are considerably higher than those typical of lines of 
sight through the disks of the Milky Way or nearby spiral 
galaxies. The only exceptions to the above are the five sys¬ 
tems at low z which have been identified with spiral disks 
and the DLA at z ~ 2.04 towards PKS 0458—020, for 
which Tg ~ 385 ± 100 K; we discuss the latter separately 
below. Note that it is not presently possible to estimate 
the spin temperatures of the remaining two L ~ L* DLAs 
at z ~ 0.437 towards PKS 1243—072 and 3C196. Finally, 
the 21cm absorption profile for PKS 1629-1-12 reported in 
this paper continues the trend (noted earlier in Chengalur 
& Kanekar (2000) and Kanekar & Chengalur (2001)) of 
low spin temperatures being found in all DLAs identified 
as L ^ L* galaxies. 

Fig. 1^ presents the data in an alternative way, 
as a scatter plot of the “true” HI 21cm equivalent 
width (i.e. (1//) / rdP) against the HI column density. 
Both of these are directly observable quantities. The 
crosses show the same quantities for Galactic HI clouds 
( polgan et al. 198^ ; Payne et al. 1982), measured using 
21cm absorption/emission studies. HI column density es¬ 
timates using observations of the Lyman-a line have been 
found to be in go od agreement with estima tes from 21cm 
emission studies ( Dickey fc Lockman 1990 ). Fig. ^ hence 
allows a direct comparison of the same observables for 
DLAs and galactic clouds. It can be clearly seen that, at 
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Fig. 4. The true equivalent width (f rdV/f) as a func¬ 
tion of the HI column density for the 24 DLAs of the 
complete sample for which we can calculate Tg values. 
Objects identihed with spiral galaxies are plotted using 
squares; the filled square shows the 2 = 0.24671 spiral 
towards PKS 1413-1-135, which may be an associated sys¬ 
tem. Other detections are shown as filled triangles, with 
non-detections shown as open triangles, with arrows. The 
crosses show measurements in the Galaxy by Colgan et 
al. (1988) and Payne et al. (1982), using 21cm emis¬ 
sion/absorption studies. 


a given HI column density, the majority of DLAs tend to 
have smaller 21cm equivalent widths than those found for 
lines of sight through the Galaxy. The only exceptions are 
DLAs that are known to be associated with spiral galax¬ 
ies (shown as squares) and the 2 ~ 2.04 DLA towards 
PKS 0458—020, (which we discuss in more detail below), 
which do lie within the range of Galactic values. The con¬ 
clusion would seem to be that the majority of DLAs con¬ 
tain larger fractions of the warm (and thus weakly ab¬ 
sorbing) phase of neutral hydrogen than typical in nearby 
spirals, as earlier pointed out by Garilli et al. (1996) and 
Kanekar & Chengalur (2001). 

It is possible that the high Tg values seen in DLAs 
might arise due to absorption originating in the outskirts 
of large disks, which would also provide a higher absorp¬ 
tion cross-section than the inner regions of these galaxies. 
Fig. I shows a plot of Tg versus impact parameter b for 
the 10 DLAs of the sample (all with 2 < 1) with optical 
identifications, as well as estimates of both Tg and im¬ 
pact parameter; absorbers identihed with spiral galaxies 
are shown as open squares. No correlation can be seen to 
exist between Tg and b. On the contrary, the largest im¬ 
pact parameter of the sample (6 = 34 kpc) is seen in a low 
Tg system, the 2 = 0.5247 DLA towards B2 0827-1-243, 
while four of the high Tg DLAs have 5 < 10 kpc. The low 


Fig. 5. The impact parameter b (in kpc) as a function 
of spin temperature for 10 DLAs at 2 < 1. The open 
squares represent objects identihed with spiral galaxies; 
the remaining detections are shown as hlled triangles, with 
non-detections shown as open triangles, with arrows. 


redshift sample clearly does not support the hypothesis 
that the high Tg seen in low redshift DLAs is because the 
absorption arises from lines of sight which pass through 
the outer regions of disk galaxies. Instead, it appears that 
the spin temperature correlates with the nature of the 
absorbing galaxy. We discuss next the 2 ^ 2.04 DLA to¬ 
wards PKS 0458—020, the only high 2 DLA with a spin 
temperature Tg < 400 K. 

The 2 = 2.03945 DLA towards PKS 0458-020 is 
as yet the highest redshift D LA with a confir med de¬ 
tection of 21cm absorption ( Wolfe et al. 1985 ). Briggs 
et al. (1989) found the VLBI 21cm absorption profile 
to be very similar to the single dish proHle and used 
this to argue that the absorber shows little structure on 
the scale of the extended radio continuum emission of 
the background source; these observations hence imply 
that the transverse size of the absorber is larger than 
8 h~^ kpc. Further, the low ionization metal lines of the 
DLA are clearly asymmetric (Prochaska & Wolfe 1997), 
which can be explained if the absorber is a rapidly rotat¬ 
ing disk (although such asymmetries are also explanable 


in other models; see, for example, Haehnelt et al. 1998 


McDonald & Miralda-Escude 1999). The spin tempera¬ 
ture of the absorber was originally estimated to be Tg = 
594 ± 297 K ( |Wolfe et al. TM^ ), the large error aris¬ 
ing due to uncertainties in htting the Lyma n-g profile 
(.Nhi = 8 ± 4 X 10^^ cm“^; Wolfe et al. 1985 ). Recently, 
however, the HI column density has been measured far 
more accurately by Pettini et al. (1994), who quote Nhi — 
4.5 ± 0.8 X 10^^ cm“^. Gombining this with the HI 21cm 
absorption profile (and taking / = 1, as confirmed by the 
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Table 4. The CNM fraction in high z DLAs. 


Name 

^ahs 

Nm 
(200 K) 
xl0^° 

Nm 

(Total) 

xl0^° 

/cnm 

PKS 0528-250 

2.8110 

<6.2 

22 ± 1 

<0.28 

TXS 2342-k342 

2.9084 

<0.7 

20 ±0.5 

<0.03 

PKS 1354-107 

2.966 

<1.3 

6.0 

<0.21 

PKS 0537-286 

2.9742 

<0.2 

2.0 

<0.09 

PKS 0336-014 

3.0619 

<0.2 

16 ± 1.3 

<0.01 

PKS 0335-122 

3.178 

<0.2 

6±1 

<0.04 

PKS 0201-kll3 

3.3875 

<1.1 

18 ±3 

<0.06 


VLBI observations) yields Tg = 385 ± 100 K, a value not 
too far from those obtained in local spirals; this is the 
first case of a DLA at z ^ 0.6 with a spin temperature 
less than 400 K. Interestingly enough, the 21cm equiva¬ 
lent width of the absorber too lies close to the range of 
Galactic values shown in Fig. |[ We note, however, that 
the metallicity of the absor ber is far lower than sola r val¬ 
ues ([Zn/H] = —1.16±0.09; Prochaska fc Wolfe 199E ), im¬ 
plying that the DLA has not undergone much chemical 
evolution. Detailed analysis of the formation of a multi- 
phase medium in proto- galaxies (^paans fc Carollo 1997 ; 


Spaans fc Norman 1997) indicate that the epoch of for¬ 


mation of a multi-phase ISM is a strong function of the 
galactic mass; galaxies with baryonic mass ~ 5 x 10^^ Mq 
form the cold phase by z ~ 3 while dwarf galaxies, with 
baryonic mass ~ 3 x 10® Mq, typically form the cold phase 
only by z ~ 1. The relatively low measured spin temper¬ 
ature in the z ~ 2.04 DLA thus makes it likely that the 
absorber is a large galaxy. The VLBI observations and the 
low spin temperature are hence both consistent with the 
DLA being a massive galaxy, with the high mass provid¬ 
ing a sufficiently high central pressure to compensate for 
the low metallicity and produce a sizeable fraction of the 
cold phase of HI. 

The evolution of spin temperature with redshift has 
been the subject of interest for some time. The original ob¬ 
servations of high spin temperatures at high redshift led to 
the suggestion that the elevation of the spin temperature 
was perhaps due to evolutionary effects. More recently, 
the det ection of low redshift DLAs (z ^ 0.2) with hig h Tg 
values (Lane et al. 1998, Chengalur fc Kanekar 199^ ) in¬ 
dicated that a simple trend did not exist in the evolution 
of Tg with redshift. However, Chengalur & Kanekar (2000) 
noted the following trend: DLAs with both low and high 
spin temperatures are found at low redshifts (z ^ 1), while 
only high Tg DLAs are found at high redshifts (z ^ 3). The 
present, larger sample of Fig. ^ shows this trend clearly, 
with the exceedingly high values of spin temperature at 
z ^ 3 indicating that these systems contain very small 
fractions of the CNM phase of HI. 


The non-detection of HI absorption in the present 
fairly deep search can be used to place stringent con¬ 
straints on the fractional content of CNM in the high z 
absorbers. To do so, we note that the allowed values of 
CNM temperatures in the Milky Way lie between 40 K 
and 200 K (Wolfire et al. 1995). We hence use 200 K (in 
Eq. (|^)) as an upper limit to the temperature of the CNM 
phase in the seven z ^ 3 DLAs to estimate upper lim¬ 
its on the CNM fraction in these systems. The original 
high resolution spectra (see Table 1) are used for this pur¬ 
pose, not the smoothed spectra used to estimate the spin 
temperature. The results are shown in Table 4 where the 
third column gives the 3tT upper limit to the HI column 
density in 200 K gas (note that this limit scales linearly 
with the assumed CNM temperature and would hence be 
even lower for lower temperature CNM). The last column 
in this table gives the upper limit to the fraction of the 
HI in CNM, /cnm; it can be seen that five of the seven 
systems at z ^ 3 have CNM fractions less than 10 %, 
with three of these having /cnm < 5 %. Clearly, most of 
the gas in the high z DLAs is warm, the typical fraction 
of CNM gas is even smaller than that in the two high 
Tg DLAs at low redshift with measured CNM fractions 
(Lane et al. 2000, Kanekar et al. 2001b) or in local dwarf 
galaxies ( Young fc Lo 1997 ). 

An independent line of evidence for high tempera¬ 
tures at high redshift comes from observations of the 
n(CH)/n(CI) ratio (Liszt 2002), which is quite different 
in the CNM and the WNM. Liszt (2002) used observa¬ 
tions of this ratio in 11 DLAs to suggest that high red¬ 
shift (z > 2.3) DLAs are dominated by warm gas, with 
at most a few per cent of the neutral HI in the CNM. 
This is in agreement with our results. Besides this, it 
has been observed that high z DLAs typically have very 
low molecular gas fractions; this can be understood in 
mode ls in which the abso rbing gas is a t high tem pera- 
ture ( Petitjean et al. 200C| ; but see also Liszt 2002 ). We 
note, in passing, that it is exceedingly interesting that the 
z = 3.0619 DLA towards PKS 0336-014 has Tg > 9240 K, 
implying that its gas is signihcantly hotter than even the 
standard WNM in the Milky Wa y (which has temp era- 
tures in the range 5000 — 8000 K; Wolfire et al. 1995 ). 

Of the 24 DLAs with Tg measurements, 16 are at z < 2 
and 8 at z > 2; we will call these the low and high z sam¬ 
ples respectively. The sample is of a sufficiently large size 
to carry out a statistical comparison to test the likeli¬ 
hood that the systems in the two samples stem from the 


same parent distribution. The Gehan test (Gehan 1965; 
[Miller 19M ) is appropriate for this comparison as it can 
be directly used in the case of limits on the measured 
quantity (in our case, the spin temperature). When ap¬ 
plied to our data, this test rules out the possibility of the 
two samples being drawn from the same distribution at 
the ^ 99% confidence level. (If one chooses to treat the 
3 ct upper limits as detections, a Kolmogorov-Smirnov rank 
1 test rules out the above hypothesis at the 95 % confi¬ 
dence level.) It is thus very unlikely that the high and 
low z DLAs are systems of similar “character”. We note 
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that 2 = 2 was used as the redshift of separation of the 
two samples as this is the approximate redshift where star 
formation is believed to peak; it is hence plausible that 
systems at 2 > 2 will show differences from those at lower 
redshifts. Since we have only 4 systems with Tg estimates 
in the redshift range 1 < 2 < 2.8 (the 2 ~ 1.4 DLAs to¬ 
wards PKS 1354-1-258 and QSO 0957-I-561A do not have 
good estimates of Tg), it is not possible to independently 
determine the threshold redshift with the present sample. 

In conclusion, the sample of damped Lyman-a systems 
which have been searched for 21cm absorption now con¬ 
sists of 31 systems; estimates of the spin temperature are 
available in 24 cases, of which 16 DLAs are at 2 < 2 and 
8 at 2 > 2. All of the low Tg, low 2 DLAs have been 
identified with large, luminous galaxies, while all DLAs 
at 2 < 1 with high spin temperature (Tg ^ 1000 K) have 
been identified either with LSBs or dwarfs. It thus ap¬ 
pears likely that the spin temperature depends mainly on 
the nature of the absorbing galaxy, as suggested earlier 
by Chengalur & Kanekar (2000). We use a recent mea¬ 
surement of the HI column density in the 2 = 2.03945 
DLA towards PKS 0458—020 to derive a relatively low Tg 
value in the absorber; this is the first high 2 DLA to show 
Tg < 400 K. The low metallicity allied with a low spin tem¬ 
perature suggest that the DLA is a massive galaxy, consis¬ 
tent with VLBI limits of its transverse size. The new deep 
GMRT 21cm spectra of this paper have also been used to 
estimate the fraction of HI in the cold phase, /cnm, in 
the absorbers at 2 ^ 3; it was found that /cnm < 0.3 in 
all seven absorbers, with /cnm < 0.1 in five of the seven 
cases. These results are in good agreement with estimates 
of .fcN M ~ few % , from observations of the n(CH)/n(CI) 
ratio ( Liszt 2002 ). Finally, we have used the above spin 
temperature estimates to test the likelihood that the 2 > 2 
and 2 < 2 DLAs are drawn from the same parent popula¬ 
tion and rule out this hypothesis at the ~ 99% confidence 
level. 
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